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" — — ^  By  means  of  easy  examples,  such  as  the  Korteweg-de  Vries,  the  Harry  Dym, 
the  sine-Gordon  equations,  and  the  Hirota  coupled  system,  it  is  shown  how 
nonlinear  partial  differential  equations  can  be  exactly  solved  by  a  direct 
algebraic  method. 

The  physical  concept,  on  which  the  method  relies,  is  one  of  generation 
and  mixing  of  the  real  exponential  solutions  of  the  underlying  linear 
equations. 

This  approach  leads  in  a  straightforward  way  to  single  solitary  waves  of 
pulse,  kink  and  cusp  shape. 

The  extension  of  the  method  towards  the  construction  of  multi-soliton 
solutions  and  the  connections  with  other  direct  methods  are  outlined.  ( 
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THE  CONSTRUCTION  OF  IMPLICIT  AND  EXPLICIT  SOLITARY  WAVE  SOLUTIONS 
OF  NONLINEAR  PARTIAL  DIFFERENTIAL  EQUATIONS 

* 

Willy  Hereman 

1.  IHTRODUCTIOH 

In  1978,  a  new  direct  method  for  generating  exact  solutions  of 
nonlinear  partial  differential  equations  (PDEs)  of  evolution  and  wave 
type  was  independently  established  by  Rosales^^  and  Korpel^2^. 

The  first  investigator  focused  on  the  mathematical  applicability 
of  a  standard  perturbation  scheme  to  construct  single  and  multi 
solitary  wave  solutions  of  many  famous  evolution  and  wave  equations. 
Korpel  looked  at  solitary  wave  formation  from  an  engineering  point  of 
view.  Applying  physical  concepts  borrowed  from  e.g.  nonlinear  optics, 
he  came  to  the  conclusion  that  a  pulse  shaped  solitary  wave  (e.g. 

O 

seek,  -type)  can  be  decomposed  into  an  obviously  convergent  infinite 
series  with  real  exponential  terms.  These  exponentials  are  nothing 
else  than  the  subsequent  harmonics  of  the  real  exponential  solution, 
characteristic  of  the  linear  dispersive  medium. 

Recently,  Hereman  et  al^3,4^  aimed  at  unifying  the  mathematical 
rigorous,  but  less  transparent  perturbation  method  and  the  heuristic 
physics/engineering  approach  toward  soliton  construction.  The 
attentive  reader  will  recognize  the  existing  isomorphisms  between 
Rosales'  iterative  scheme  and  our  recursive  system,  between  the 
solution  techniques  and  the  summation  procedures  applied  to  the 
resulting  infinite  series  expansions. 

In  this  paper  we  illustrate  by  means  of  rather  easy  examples  how 
the  physical  interpretation  deepens  the  understanding  of  solitary  wave 
formation.  Lack  of  space  does  not  permit  us  to  elaborate  on  the 
complementary  mathematical  details,  which  can  be  learned  from  earlier 
work^1,4,S^  on  the  subjeot. 
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Advocating  a  pedagogical  approach,  ve  first  construct  the 
well-known  hyperbolic  secant  squared  solution  of  the  Korteweg-de  Vries 
equation  (KdV),  working  our  way  up  to  the  derivation  of  the 
multi-soliton  solution  of  the  sine-Gordon  equation  (sG)  as  a  finale. 

In  between,  we  solve  the  Harry  Dym  (HD),  the  coupled  Korteweg-de  Vries 
CoEdV)  and  the  sine-Gordon  equation,  for  which  the  construction  of 
single  solitary  waves  is  of  intermediate  level  of  difficulty. 

Wherever  appropriate  we  indicate  connections  with  various  other 
direct  methods  (Trace  method''  ,  Sirota's  bilinear  formalism  J  and  its 
clones^6, 9\  Fredholm  determinant  method^19 ’ ,  direct 
linearization^12^,  direct  integration^13^,  etc.).  In  this  manner  we 
outline  the  general  framework  in  which  this  contribution  fits. 

2.  THE  KORTEVEG-DE  VRIES  EQUATIOH 
The  celebrated  KdV ( 14-16 ^ , 

Ut  +  auu^  +  U3x  -  0  ,  a  €  R  ,  (1) 

3  3 

wherein  subscripts  indicate  partial  derivatives  (e.g.  u3x  -  d  n/d  x), 
describes  shallow  water  waves,  ion-acoustic  waves  in  plasmas,  the 
dynamics  of  a  nonlinear  lattice,  eto. 

Using  (1)  as  a  paradigm,  we  sketch  our  direct  algebraic  method  in 
ten  steps,  confining  ourselves  to  a  single  solitary  wave. 

(i)  Searching  for  a  stationary  solution,  we  introduce  a  moving  frame 
of  reference  by  assuming  u(x,t)  -  ♦Cf ) ,  where  f  -  x  -  vt. 

The  constant  v  refers  to  the  anticipated  travelling  wave 
velocity.  Doing  so,  PDE  (1)  can  be  replaced  by  the  ODE 


(ii)  For  mathematical  convenience,  we  next  integrate  (2),  yielding 

-v*  +  (a/2)*2  +  +  cC  -  0  ,  (3) 

where  cC  is  some  integration  constant. 

(ill)  Ve  will  allow  travelling  wave  solutions  to  have  a  constant  term 

a 

o.  Upon  substitution  of  *  -  o  +  *  into  (3),  we  get 

(aC-v)*  +  (a/2)*2  +  +  c(-V+(aO/2)+C)  -  0  .  (4) 

r  o') 

(iv)  Adhering  to  the  intuitively  simple  model  of  solitary  wave 
formation,  we  solve  the  linearized  version  of  (4),  (i.e.  ignore 
the  second  term)  for  real  exponential  solutions  in  the  form 

♦  (f)  -  exp(k(v)f).  Obviously  the  dispersion  law, 
k2(v)  -  v  -  ao  -  2C  -  v,  follows  provided  c  -  (2/a)(v-C). 
Observe  that  both  decaying  and  rising  exponentials  are  possible. 
For  the  remainder  we  will  work  with  the  decaying  solution, 
denoted  by 

g(f)  -  exp(-kf),  k  -  V2C-v  .  (5) 

(v)  Appropriate  scaling  by  ♦  -  (2/a)(2C-v)*  -  (2k2/a)*,  simplifies 


(4)  into 


-♦  +  *2  +  a/k2)*2f  -  o  .  (6) 

(o') 

(vi)  According  to  Korpel's  model  ,  the  quadratic  nonlinearity  in 
(6)  will  square  the  linear  solution  g.  Subsequently,  g  and 

O 

g  now  being  present,  in  the  next  step  the  nonlinearity 
3  4 

generates  g  and  g  .  Proceeding  in  this  way,  any  integer 
power  of  g  will  be  created.  In  the  optics  terminology, 
angn  -  aQexp(nkx-nut ) ,  with  wave  number  k  and  angular 
frequency  <j  -  kv,  represents  the  ntlx  harmonio  wave  with 
amplitude  an-  For  n  -  1,  we  refer  to  a1g  as  the 
fundamental  (wave).  This  heuristic  principle  gave  Impetus  to 
the  search  for  an  exact  solution  (to  the  nonlinear  equation)  in 
the  form 


♦  -  l  aaAf)  • 

n-1 

(vii)  To  determine  the  coefficients  aQ,  we  substitute  (7)  into  (6) 

(3  4) 

and  apply  Cauchy's  rule  *  to  recollect  equal  powers  in  g. 
This  results  in  an  Infinite  hierarchy  of  equations: 

n-1 

(n2-l)an  +  J  a£an-e  "  °’  n  *  2  •  (8) 

l-l 

starting  with  an  arbitrary  (positive)  constant  a^  To  the 
physicist,  the  recursion  relation  (8)  for  the  amplitudes  merely 
describes  the  mechanism  of  energy  transfer  into  successive 
modes. 


-4- 


(vill)  Solving  the  system  (8)  recursively,  ve  recognize  that 


an  -  6n(-l)n+1an,  a  -  ax/6  >  0,  n  €  »  .  (9) 

In  an  earlier  paper  we  discussed  various  techniques  to  solve 
systems  similar  to  (8),  how  to  recognize  the  general  solution 
and  how  to  verify  it. 

(ix)  The  last  but  one  step  aims  at  the  summation  of  the  resulting 
infinite  series 

Oft 

♦  -  6  J  C-l)n+1n(ag)n  ,  CIO) 

n-1 

which  converges  only  for  sufficiently  large  f,  into  its  closed 
form 

♦  -  Oag/Cl+ag)2  ,  (11) 

which  is  valid  everywhere . 

(x)  Finally,  writing  (11)  in  the  original  variables  x  and  t, 
gives 

u(x,t)  -  (2/a)(v-C)  +  (3/a)(2C-v)sech2{^V2C-v(x-vt)  +  6}  .  (12) 

The  single  solitary  wave  solution  (12)  depends  on  three  arbitrary, 
unrelated  constants  v,C  and  6  -  g£n(l/a).  Observe  that  for  the 
particular  choices  C  -  v  >  0  and  C  -  v/4  <0,  we  obtain  the 
familiar  pulse-type  solution, 


(13) 


u(x,t)  -  (3v/a)sech2{j£>/v(x-vt)  +  6}  , 
and  well-type  solution 

u(x,t)  -  (3v/2a)tatnh2{^V-v/2(x-vt)  +  5}  .  (14) 

Of  course,  finding  the  solitary  wave  form  (12)  is  neither  original  nor 

a  great  achievement;  it  has  been  obtained  by  half  a  dozen  different 
(14-17) 

methods  ' ,  but  this  method  clearly  reveals  the  physical  mechanism 

behind  solitary  wave  formation. 


3.  THE  HARRY  DYM  EQUATION  J 

The  prototype  of  equation  for  so-called  cusp  (or  spiky)  solitary  } 

wave  solutions  is  the  Harry  Dym  equation^ 10-21  \  jj 


ut  +  (l-u)Su3x  -  0  ,  (15) 

(19) 

or  any  equivalent  f ormv  ' ,  which  occurs  in  connection  with  the 
classical  string  problem. 

The  presence  of  a  cubic  nonlinearity  in  the  coefficient  of  the 

dispersive  term  u3x  drastically  changes  the  nature  of  the  most 

(20  21) 

elementary  particular  solution  ,  i.e. 


(16) 


u(x,t)  -  sech2{|v*v[x-vt+6(x,t)] }  , 


where  the  time  and  space  dependent  phase  6  is  governed  by  the 


transcendental  equation 


6(x, t )  -  c  +  (2/V!v)tanh{^/v[x-vt+6(x, t) }  , 


(17) 


where  c  is  an  arbitrary  constant.  This  implicit  solution  was 
obtained  through  the  inverse  scattering  technique  (1ST). 

Straightforward  application  of  our  method  to  (15),  which 
accidentally  has  the  same  linearized  version  as  the  KdV  (1),  does  not 
lead  to  any  solution.  Obviously,  we  must  relax  our  approach  by 
introducing  the  new  variable. 


f(x,t)  -  k[x-vt+6(x,t)3  , 


(18) 


and  then  search  for  solutions  u(x,t)  -  F(f),  with  k6(x,t)  -  G(f), 
where  both  F  and  G  remain  to  be  determined. 

Through  the  operator  relations 


a.  -kv  d.  a.  k  d. 
at  "  df  ’  ax  "  1-Gf  df 


(19) 


(15)  is  transformed  into 


-kvFf  +  k3(l-F)3(l-Gf)~4{(FfG3f+3F2fG2f)(l-Gf) 


♦  F„„(l-G„)2  +  3F*G2  >  -  0  . 


(20) 


>  «.•  I 


Clearly,  Gf  -  P  is  the  right  choice  to  simplify  (20)  to 


(-v/k2)(l-F)Ff  +  [4FfF2f(l-F)+F3f(l-F)2+3F^]  -  0 


(21) 


After  division  throughout  by  (1-F)  ,  followed  by  a  first  integration 
(introducing  constant  c.^) ,  a  subsequent  multiplication  by  Ff  and 
another  integration  (constant  c2),  we  arive  a<b 


F2  -  (l-F)[(v/k2-c2)+(c2-c1)F+c1F2]  . 


(22) 


Seeking  for  a  solitary  wave,  we  add  the  boundary  conditions 
F,Ff,F2f  -»  0  as  |f  |  Therefore,  set  o1  -  c2  -  v/fc2-  and 

readily  integrate  (22),  yielding 


F(f )  -  sech2{>/v/2k)(f+c3)}  , 


(23) 


where  c3  is  the  final  integration  constant.  With  the  definition  (18) 
of  f,  we  thus  obtain  solution  (16).  Regarding  the  choice  Gf  -  F, 

G  -  k6;  after  integration  of  (23)  we  retrieve  (17),  wherein 
c  -  (c3+c4)/k  and  c3  is  absorbed  in  5.  Other  relevant  solutions  of 
(22),  as  listed  e.g.  by  Drazlnv  J  are  presently  under 
investigation^19^ . 

Although  we  decided  in  favor  of  direct  integration  of  (21), 
application  of  our  direot  method  would  effortless  have  led  to 
(16)-(17).  Hence,  in  conclusion,  a  slight  generalization  of  our  method 
broadens  the  class  of  retrievable  solutions  to  implicit  ones. 


'  v.v-y:o:w 


> 


ww. 


4.  THE  COUPLED  KORTBVEG-DE  VRIES  EQUATIONS 

To  extract  further  information  about  the  applicability  of  our 


technique  let  us  investigate  how  it  would  fare  on  a  famous  coupled 
system^23-29^  r 


u^  -  “(6uux+u3s;)  -  2|3wwx  -  0  ,  (24) 

W*  +  3uw_  +  W_  -  0,  a,p  €  R  .  (25) 

u  X  OX 

These  equations  describe  the  interaction  of  two  long  waves  with 

(24  25) 

different  dispersion  laws.  The  coupled  system  is  often  referred 

to  as  the  coupled  KdV  equations,  as  for  w  -  0  it  reduces  to  the 

KdV  in  u.  Sometimes  (24)-(25)  is  quoted^26-28^  as  the  Hlrota-Satsuma 

system  after  the  two  investigators  that  first  solved  it  using  a  quite 

ingenious  bilinear  formalism^  *  . 

(29) 

In  a  forthcoming  paperv  J  we  will  prove  that  if  u  is  of 
travelling  wave  form,  say  u(x,t)  -  *(f),  with  f  -  x  -  vt,  then  w 
exhibits  the  same  form,  hence  w(x,t)  -  *(£).  Thus  we  must  solve 

V*  +  3a*2  +  a*-.  +  0*2  -  0  ,  (26) 

-  +  3**f  +  -  0  ,  (27) 

ignoring  integration  constants  here.  Substitution  of  g(f)  -  exp(-kf) 

into  the  linear  parts  of  (26)-(27),  leads  to  two  dispersion  laws, 

2  2 

v  -  -ak  and  v  -  k  .  According  to  our  philosophy,  the  nonlinear 


-9- 


solutions  <t>  and  ^  can  only  be  built  up  from  the  same  real 

(29) 

exponential  if  a  -  -1.  A  profound  study  reveals  that  either  one 
of  the  dispersion  laws  leads  to  the  same  exact  solutions.  Hence,  let 
us  proceed  with  k  -  Vv,  v  >  0. 

Upon  substitution  of  the  scaled  series  representations, 


*  "  I  1  angn(f),  + 

n-1 

into  (26)-(27),  we  arrive  at 


l  bng a(f)  , 

^TpT  nfi  n 


(28) 


(l+an2)an  +  Y  ^a&£&n-e+ehehn-e^  ”  °’  n  ^  2  . 


n(n2-l)bn  +  J  ehe&n-e  “  0  *  n  1  2  , 


(29) 


(30) 


together  with  (1+cOa^  -  0,  b^  arbitrary,  and  where  e  -  +1  (e  -  -1) 
if  p  >  0  CP  <  0). 

Two  interesting  cases  pop  up: 

(i)  For  a  jt  -1 ,  so  a1  -  0,  we  straightforwardly  obtain 


a2n  "  24(-l)n+1nAn,  b2a  -  0.  n  €  N\<0>  , 

(31) 

a2n+l  -  °-  b2a+l  -  '-“V'  “  e  "  ■ 

(32) 

f4 

■a 

■A 

L 

S 

8 


requiring  A  -  -eb2/24(4a+l)  >  0,  hence  p  and  4a  +  1  must 
have  opposite  signs.  With  (31)-(32)  the  solutions  (28)  then 


become 


♦  -  av  y  (-Dn+in<AgV  -  -ev-^: , , 

Cl-Ag2)2 


(33) 


+  -  -z—  )  c-dVaV”*1  -  — - hr  • 

v^TFT  n:0  1  V3JFT  (i+Ag2) 


Returning  to  the  original  variables,  we  obtain 


(34) 


u(x,t)  -  2v  sech2[Vv(x-vt)+6]  , 


w(z.t) 


i/p  sech[Vv(x-vt)+6]  , 


(35) 


(36) 


With  6  -  |24(4a+l)/b2 | .  Note  that  the  cKdV  indeed  remains 

invariant  for  reversing  the  sign  of  w.  For  a  -  -1/4,  clearly 
bQ  «  0,  Vn  c  N;  so  w  -  0  and  the  cKdV  reduces  to  the  KdV  (1) 
with  a  -  6.  Replacing  t  by  t/4,  v  by  4v,  from  (13)  we  get 
(35).  We  should  remark  that  (36)  is  a  special  case  of  a  solution 
obtained  by  the  dressing  operator  technique^27^ .  Furthermore,  for 

(  OR") 

a  -  1/2  the  cKdV  are  known  to  be  completely  integrablev  ,  it 
has  the  Painleve  property^26^  and  a  N-soliton  solution^24^ . 


■rare* 


(ii)  For  a  -  -1,  apparently  a^  is  arbitrary  and  so  is  b^. 

Recursively,  one  can  calculate  all  an  and  bn  in  tbe  hope  to 
obtain  the  general  closed  form,  which  we  do  not  know  yet. 


2  2 

However,  for  b1  -  a1/2 

and  p 

>  0 

(so  e  -  1) 

we  obtained 

a  -  12n(-l)n+1an,  b2  - 
n  n 

4/z- 

with 

a  -  a1/12. 

Substitution 

into  (28),  leads  to 

u(x,t)  -  v  sech2[|Vv(x-vt)+6]  (37) 

v(x,t)  -  (3/V6  \(5  ( )  u(x,t)  -  (3v/V6  |/3  | )  sech2[^Vv(x-vt)+5] ,  (38) 

with  6  -  ^  «#i(12/a1).  Observe  that  for  w  -  3u>/sp  both 
equations  (24)-(25)  become  identical  to  the  KdV  (1)  with  a  -  3, 
(13)  then  being  the  same  as  (37). 

5.  THE  SUBTK-GORDOH  EQUATION 

The  sG  equation,  in  light  cone  coordinates, 

u^  -  sin  u,  (39) 

describes  the  propagation  of  crystal  dislocation,  superconductivity  in 
a  Josephson  Junction,  ultrashort  optical  pulse  propagation  in  a 
resonant  medium,  etc. 16,22\  For  the  mathematician,  (39)  is  long 
known  in  the  differential  geometry  of  surfaces  of  constant  negative 
curvature^ 15 ’ 16 ^ . 


At  the  cost  of  dealing  with  one  equation  with  a  transcendental  ( ! ) 
nonlinearity,  we  rather  transform  the  sG  into  a  nonlinear  coupled 
system  with  strictly  polynomial  terms: 

o  .  (40) 

2*  +  #2  +  $2  -  0  ,  (41) 

where  *  -  u  ,  +  -  (cos  u)  -  1  . 

To  construct  a  single  solitary  wave,  we  proceed  as  in  section  4, 
focusing  on  steady  solutions  $(f)  -  $(x,t),  -  *(x,t)  with 

f  -  x-vt.  Expanding  the  scaled  functions  as 

00  09 

♦  -  (l/v^v)  l  a ngn(f),  *  -  l  bngn(J:)  *  C42) 

n-1  n-1 

with  g(j;)  -  exp(-kf),  we  obtain 

n-1 

(n2-l)an  -  J  *e\-e  m  0  •  n  2  2  .  (43) 

«-l 

n-1 

2bn  +  1  (b<V<+Kn_<)ft<an-P  "  0  ■  11  *  2  •  (44) 

«-l 

2 

where  we  used  the  dispersion  law  k  -  -1/v,  v  >  0,  to  simplify. 
Iterative  calculation  suggests 


^^Bi 

I  *  ’  i 

a2a  -  0.  b2n  -  8(-l)nna2n.  n  e  N\{0}  . 
a2a+l  "  4^-1^  a  »  ^n+l  “  n  €  N  * 


where  a  -  a1/4  >  0.  Upon  substitution  into  (42),  we  find 


n,  „,.2n+l  4vcv  a« 


H 

(45) 

H 

(46) 

liB** 

■l 

*  -  (4/V=v)  l  (-l)n(ag) 
n-0 


l+(ag)‘ 


(47) 


♦  -  -8  J  (-l)n+1n(ag)2n  -  .zB&SgL  , 
nil  [l+(ag)2]2 


In  the  variables  x  and  t,  we  thus  get 

(cos  u(x,t))  -  1  -  1  -.2  se<5h2[(l/>^v)(x-vt)+5]  , 

u(x,t)  -  i  — —  f  sech[(l/-/:v)(x-vt)+6]dx 
V^v  J 

-  ±4  arctan{exp[(l/'/:v)(x-vt)+6] }  , 


(48) 


(49) 


(50) 


with  5  -  «n(4/a1).  This  is  the  well-known  kink-type  solution  of  the 


sine-Gordon  equation 


(11,15,16) 


6.  H-SOLITOH  SOLUTIONS 


The  most  effeotive  techniques  to  construct  N-soliton  solutions 
are  inverse  scattering^15-18, 22^  and  Hirota's  method^7, 15 ’ 16^ ,  the 
latter  being  olosely  linked^1 ’ 28 ’ to  all  other  (often  iterative) 
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procedures  listed  in  the  introduction.  Neither  of  these  methods  is 
very  transparent  in  explaining  why  N-solitons  are  the  way  they  are. 

(  A} 

In  the  spirit  of  an  earlier  paperv  ,  where  we  showed  how  the 
N-soliton  solution  to  the  KdV  could  be  built  up  from  N  real 
exponentials,  we  construct  the  N-soliton  solution  of  the  sG  (39), 
using  (40)  and  (41).  Emphasis  is  again  on  physical  relevance,  more 
than  on  mathematical  rigor. 

Motivated  by  the  fact^  that  multi-soliton  solutions  evolve  from 
the  mixing  of  real  exponential  solutions  of  the  underlying  linear 
equation,  we  substitute 

N  N 

“  I  Oig1(x,t)  -  J  c^  exp(k1x-toit)  ,  (51) 

i-1  i-1 

into  the  linear  part  of  (40),  implying  -  -1/k^  (i  -  1,2, •••,N). 

The  constants  ci(ki)  will  be  fixed  later.  Observe  that  the  starting 

(2) 

term  in  the  expansion  of  tf,  say  ¥  ,  must  be  of  the  form 


so  that  -Z'i 


balances  the  term 


we  find 


in  (41).  Thus,  with  the  dispersion  law  -  -1/1^ 


1  1  °i°i  l  V"! 

diJ  _  ~2  ci0JuluJ  '  “2  “  2  °i°J  k±+kj  * 


(54) 


Note  that  any  term  in  the  expansion  of  4>  (*,  respectively)  will  only 
have  an  odd  (even,  respectively)  number  of  g's,  which  is  in  agreement 
with  (47)  and  (48). 


The  analogue  of  the  first  term  in  (43),  i.e.  (n  -l)an  ,  will 
result  from  the  action  of  the  linear  operator 


L*  “  fiat  - 


(55) 


on  the  (2n+l)th  term  in  the  expansion  of  $,  namely 


,(2n+l) 


N  N 


I  l  ”•  l  *1*J  ”•  •  n  £  NX<01 


(56) 


,i-l  J-l^  B-j 
2n+l  summations 


The  analogue  to  the  second  term  in  (43)  will,  in  its  most  symmetric 
form,  look  like 


5  1 


,(2«+l) 


(k^.k, ,  •  •  •  ,kQ)  *(2n-2«)  (k  ,k  .•••.ks) 


2£+l  arguments 


2(n-  e)  arguments 


.  ^(2n-2^  )  /t,  y  y  \  - (2£+l)  /y  y  y  \ 

+  v  (k^.k, , • • • ,xg)  (kB,kn, * • * ,kg) 


(57) 


2(n-«)  arguments 
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The  analogue  to  (44)  reads 


*(2n)  -  l  l 


N  N  N 

I  d±J • • *  r  ^i^4  *^r 

i-1  j-l  r-1 


2n  summations 


-  l  *(an-20^o^".y 

2£  arguments  2(n-£)  arguments 


n-1 


+  l  ♦Ja<+l)rt1,k,.....k0)  v ...,kE).  n  2  1. 


£-0 


2£+l  arguments 


2(n-0-l  arguments 


allows  to  subsequently  determine  the  coefficients  d^  ><r  . 

To  make  this  less  obscure,  let  us  give  an  example,  and 

(2) 

#  being  computed,  we  equate 


N  1  N 


W  3  '  l  \  1  C-(ui+uJ+uk)(VkJ+kk)_1I°iJk  *i*j*k 

i-1  J-l  k-l 


N  N  N 

’III  (VV(VV(Vkl>(°ljk/klkjV*l*J«! 

i-1  4-1  k-l 


to 


-17- 


(58) 


(59) 


+  +C2)(ki,kJ)0(1)(3tk)3 

H  N  K 

"  “5  I  1  1  (ki+kJ)(kj+ltk)(ltk+]ci)(cicJck/ltikjkk)gigjgk  *  (60) 

1-1  J-l  k-1 

therefore  c^^  “  -1/ t4(k1+kj )(kj+kfc) ]  if  we  set  -  1.  After  some 
lengthy  calculations,  partly  carried  out  with  MACSYMA  (i.e.  a  large 
scale  computer  program  that  performs  algebraic  manipulations),  we 
obtain 


I 


K* 


...s  -  C-Dn4  n[(ki+k  )(k J+kk)--.(kp+kB)]_i.  n  e  n. 


dl J •  •  •  r  “  ("1)n  lgl  2nC‘*>i+«j+“-+wr)Uki+kJ)CkJ+kk)*--(kq+kr)]"1 


n  €  N\{0>.  (61) 


The  final  objective  is  then  to  write 


,(2n+l) 


N  N  N  N 

l  l  l  Tsgfe*?  + 

L-l  1-1  J-l  k-1  J  J  * 


N  N 

I  1 

i-1  J-l 


.  X  (_ nD/hn  gjgj  *•*  Ss _ 


.  .  (62) 


and  a  similar  expression  for  in  their  closed  forms.  Various 

authors^ 11 .31 )  ^ave  shown  that  this  can  be  done  by  introducing 
the  N  x  N  matrices  I  (unity)  and  B,  with  elements 


The  N-soliton  solution  to  the  sine-Gordon  equations  (40)-(41)  is  then 
found  to  be 

4>(x,t)  -  4[Tr(arctan  B)]^  , 

*(x,t)  -  -2t«ft(det(I+B2))]xt  . 

while,  regarding  <*>  -  ux  , 

u(x,t)  -  ±4  Tr(arctan  B)  -  ±2  Tr  [fzxf]  ]  (66) 

satisfies  (39). 
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